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Application of modified hat functions 
for solving nonlinear quadratic integral 
equations 


F. Mirzaee* and E. Hadadiyan 


Abstract 


A numerical method to solve nonlinear quadratic integral equations (QIE) 
is presented in this work. The method is based upon modification of hat func- 
tions (MHFs) and their operational matrices. By using this approach and 
the collocation points, solving the nonlinear QIE reduces to solve a nonlin- 
ear system of algebraic equations. The proposed method does not need any 
integration for obtaining the constant coefficients. Hence, it can be applied 
in a simple and fast technique. Convergence analysis and associated theo- 
rems are considered. Some numerical examples illustrate the accuracy and 
computational efficiency of the proposed method. 


Keywords: Modification of hat functions; Nonlinear quadratic integral 
equation; Vector forms; Operational matrix; Error analysis. 


1 Introduction 


Over the last years, the integral equations have been used increasingly in 
different areas of applied science. This tendency could be explained by the 
deduction of knowledge models which describe real physical phenomena. For 
details, we refer to [1, 2, 4-9, 12-18, 23, 25, 26]. In particular, quadratic 
integral equations (QIEs) have many useful applications in the real world. 
For example, QIEs are often applicable in the theory of radiative transfer, 
the kinetic theory of gases, the theory of neutron transport, the queuing 
theory and the traffic theory. The QIEs can be very often encountered in 
many applications. The quadratic integral equations have been studied in 
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several papers and monographs [1, 2, 4-26, 29, 30]. In this paper, we study 
the numerical solution of a QIE: 


se) = ae) + (f° entity. repay) (f° ba(en)Ualu fonay) 
(1) 
where x € D = [0,1], g(x) € C3(D), Ui (2, f(x)), Ue(a2, f(x)) € C3(D x R) 
and ki(a,y), ko(x,y) € C3(D x D) are known functions, f(x) € C3(D) is 
an unknown function and we will obtain the approximate solution in the 
truncated MHFs series form 


so that fi; ¢ = 0,1,...,m, are the unknown MHF’ coefficients and h;(x); 
i=0,1,...,m, are the MHFs. 

To mention some recent works on QIEs, see e.g., [1,8,12, 14,15, 25,26] and 
for some applications we refer readers to [18,23]. Existence, uniqueness and 
some other properties of the solution to these problems were established in 
[33]. It should be recalled that nonlinear QIEs have been treated extensively 
with the measure of noncompactness and a fixed point theorem of Darbo 
type. This approach seems to be too restrictive. Furthermore, in most of the 
above investigations, some additional assumptions in terms of the measure 
of noncompactness were imposed on g(x). 

The plan for this paper is as follows: In Section 2, we describe MHFs 
and their properties. In Section 3, we will apply these sets of MHFs for 
approximating the solution of QIEs. In Section 4, theorems are proved for 
convergence analysis. Numerical results are given in Section 5 to illustrate 
the efficiency and the accuracy of our algorithms. Finally, Section 6 concludes 
the paper. 


2 Modification of hat functions 


The purpose of this section is to collect a number of definitions and lemmas 
concerning MHFs. we first construct the set of MHFs. 


An (m+ 1)-set of MHFs consists of (m+ 1) functions which are defined 
over district D as follows [3,29] 


az (x —h)(a— 2h) O< a < 2h, 
ho(x) = 


0 otherwise, 


ifiisodd andl <i<m-—l, 
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s(x — (i- 1)h)(e@- (6+ 1h) GG -DA <a t+ Da, 
0 otherwise, 


if zis even and2<i<m-2, 


spa (@ — (i -1)h)(@—- (4 -2)h) (6 -2)K <M a<ih, 
hi(z) = 4 spx(a@— (6+ Ih)(x@— (6 +.2)h) th x < (i +2)h, 
0 otherwise, 
and 
syz(u —(1—h))(w@- (1-2h)) 1-2h<aK<l, 
Pm (@) = 


0 otherwise, 
where m > 2 is an even integer and h = 4. It is obvious that 
li=j, 
hi(gh) = (2) 
O1F3), 
0 7 is even and |i — j| > 3, 
0 4 is odd and ji — j| > 2, 


and 


Let us write the MHFs vector H() as follows 


H(x) = [ho(z), hi(a),...,Am(a)|"; 2 € D. (4) 
An arbitrary function f(x) defined over D can be expanded by the MHFs 
as 
f(x) ~ FT H(x) = H7(2)F, 
where 
F= [fostinertn tel 5 
and 


fi = f(ih); i=0,1,...,m. 


Similarly an arbitrary function of two variables, k(x, y) on district D x D 
may be approximated with respect to MHFs such as 
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k(a,y) ~ H" (2) KH(y), 


where H(x) and H(y) are MHFs vector of dimension (m+ 1) and K is the 
(m+ 1) x (m+ 1) MHFs coefficients matrix. 


According to (2) and expanding te hi(y)dy, i = 0,1,...,m by MHFs, 
integration of vector H(x) defined in (4) can be expressed as 


iu  H(y)dy © PH(c), (5) 


where P is the (m +1) x (m+1) matrix as follows 


ee ae ae ae eee ae ae 
0 8 16 16 161616... 16 16 
0-149 888...8 8 

nh | 0 0 8 16 161616... 16 16 
P=75)/000-1498...8 8 
0000000... 8 16 
0000000...-14 


By considering (2), (3) and expanding entries of H(x)H* (x) by MHFs, 
we obtain 


H(a)H"™ (x) ~ diag(H™ («)), 
so we have o 
H(«)H" («)F ~ FH(a), (6) 
where F be an (m+1)-vector and F is an (m+1) x (m+1) diagonal matrix. 
Also, if A is an (m +1) x (m+ 1)-matrix, we have 


H™ (x) AH(x) ~ H™(2x)A, (7) 


where A is an (m + 1)-vector with elements equal to the diagonal entries of 
matrix A. 


3 Basic idea 


In this section, we will provide the basic idea. This idea includes of 
approximating the solution of nonlinear quadratic integral equations (1). To 
solve this equation, we first consider the approximations 
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wi(x) = Ui (2, f(x)) 
=U, (a, g(a) + (fp ki(a, y)wi(y)dy) (Jo ko (a, y)we(y)dy)) , 
w(x) = Ua(a, f(x)) 


= U2 (x, g(x) + (fo k(x, y)wi(y)dy) (fo a(x, y)wa(y)dy)) - 


(8) 


We approximate function w1(x), wo(x), ki(a,y) and ke(x,y) by MHFs, 
wi(«) ~ WP A(x) = HT (2), 
w(x) ~ Wy A(x) = H7(2)We, 
ki(v,y) ~ H"(«)K1H(y) = H"(y)K17H(a), 
ko(v,y) ~ H" (x) K2H(y) = H" (y)K27H(a), 


where Wi, W2, K1 and K2 are MHFs coefficients of wi(x), w2(x), ki (x, y) 
and ko(x, y), respectively. Substituting (9) in (8), we get 


HT (e)W, = Ui(2,9(0)+ ( fp H7(e) 1H (y)H? (y)Widy) 
x (fg) H? (x) K2H(y)H" (y)Wady) ), 

HT ()W, = Us(<, g(@)+ (fy HT(e)K1H(y) HT (y)Widy) 
x ( fy H? (a) K2H(y)H? (y)Wody) ). 


Using (5) and (6), yields 


HT («)W, ~ 0 (2, g(x) + (#7 (e)K1W,PH(e)) (H7(2)K2W2PH(x))), 


HT (x)We ~ Ur (x, 9(@) i (#7 (©) K1W,PH(c)) (7 («)K2W2PH(c)) ) 


where W; = diag(W;), i= 1,2, are an (m +1) x (m+ 1) diagonal matrices. 
From (7), we have 


ATW & v1 ( ree (17 (ey Ku07,P) (w"(e)K2IP) ), i” 


HT (2)Wo ~ Ua ( g(x) + ("e107 P) (w"e)K2WaP) ). 
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eee 


where KiW;P, i = 1,2, be an (m+ 1)-vector with elements equal to the 


diagonal entries of matrix Kiw;P. We can rewrite Kiw,P, 7 = 1,2, as 
follows 


KiW,P = AiW;, i=1,2, (11) 
where 
Aipg = KipgPap, p,¢=9,1,...,m 


Substituting (11) into (10) and replacing ~ by =, we obtain 


H? (x)W, =U, (2, g(x) + H7 (x) A1W, H? (x) A2W3) , 


H? (x)W2 = U2 (2, g(x) + HT (2) AW, H7 (x) A2W2) . 
Now, using Newton-Cotes nodes as 


21 —1 
2(m+1)’ 


LS 


then 


H? (x;,)W, = Uj (24, 9(xi) + H? (x;)A1W, H7 (2;) A2W2) , 


H? (x;)W2 = Up (24, g(x) + H? (x;)A1W H7 (2;) A2W2) . 


We have a system of (m+1)? nonlinear equations and (m+ 1)? unknowns. 
After solving the above nonlinear system, we can find W; and W2 and then 


f(x) ~ fm(x) = g(x) + H7 (x) A1W,H™ (x) A2W2. 


4 Convergence analysis 


In this section, for confirming the accuracy of the proposed scheme in 
the previous section analytically, we provide an upper bound for difference 
between the exact solution of (1) and our approximated solution. We show 
that the MHFs method, is convergent of order O(h?). We define 


IIgl| = sup |g(x)]. 
xE€D 


Theorem 1. Suppose x; = ih, i=0,1,.. g(x) 
the MHFs expansions of g(x) that defined as ate = 
assume that em = ||g — gm|| where x € D, then 


o* (D) and gm(a) be 


E 
9 9 (a1) hile). Also, 


Cm = O(R*): 


Application of modified hat functions for solving nonlinear ... 71 


Proof. According to definition of MHFs, the g(a) is the quadratic poly- 


nomial interpolation on [x;-2,x;). Therefore, for the interpolation error on 
[x;-2,%;), we have [31] 


pice, 
Ci jn( 2) = 9(£) — Gm (av) = zd I&4) Il (a—ay); t= 2,4,...,m 


6 dx 4. 
uf =i-2 


where «, &; € [x;_2, ai). Let u(x) = [[,_,_2(2 — air), so 


1 9 9(Ei) 
lleim||l == sup v(a)} , +=2,4,...,m, 
6 LE [xj—2,24) os 
° d9(E) 
1 VS 5 
les mll Ss in sup | | ( )| ’ a 2,4, »™m. 
6 x€[xj~2,2;) dx* 
On the other hand, we have 
Cm = sup |g(x) — gm(x)| = max sup —|g(x) — gm(2)| 
2e€D $= 2A es ™ xE[x;_2,%;) 
= max lem 
1=2,4,..., 
So ; p 
em <=. max ss (6%) lv(x)| - 
6 t= 2A, m x€[v;_2,2%) dx 


Since |v(x)| < supze; 


of ivan 2 (x = it) 


Tice — 241) 


is obtained at x = (i-1—- Y3)h, we have 


i es and the maximum value 


lu(z)| S 


ee Va © [aj_-2, 0). 
Therefore, it is not difficult to verify that 
3 
em < < is z= (12) 
So 


em = O(h?). 


This completes the proof. 


Theorem 2. Let x; = y; = ih, i=0,1,...,m, k(x,y) € C3(D x D) and 
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i=0 j=0 
be the MHFs expansions of k(x,y). Then 
em = O(h*), 


where €m = \|k — km || and (a,y) € D x D. 


Proof. km(«,y) is the quadratic polynomial interpolation of k(x, y) on Q4; = 
[v;-1, 2) X [yj-1,y;). Therefore for the interpolation error on 9;;, we have 
[27] 

€ij,m(2,Y) k(a, y) ~ km (a, y) 


_ 1 0%K(&,y) ieG=eee 1 OPk(a,n;) Il ow 


6 Ox? 4. 6 dy? a 
i! =1-—2 ji =j-2 
1 O8k(E,m) 7 ! 
~36 on307 LI @-) TT ww), 
i! =i-2 j!=j-2 


where a9 _ 2,4, ree ™M, Dye, € [vi-1, Zi) and YN; € reeeatrpe So we 
have 


-2 Pk(E, y) _ Bk(a,n;) 1 O°K(E, 5) 
\|€:3,m| a 6 Gapees Ax3 v(a) T dy? u(y) 6 d230y3 vu(a)u(y) ) 
or 
1 Pky) k(x, ny) 
Me aes ie WRU de poceies oee E 
loom <6 ame {| fee] joe |S] tata) 
1 OPk( 1) 
6 Bap lv(x)| Ju(w)l ep. 


where i,j = 2,4,...,m, v(x) = TI —;_2(#—-24) and u(y) = Th =j-2(y-yy") 
On the other hand, we have 


Em = sup |k(x, y) — km(x, y)| 
(x,y)EDxD 
= max —sup__|k(z,y) —km(2,y)|=_ max — |leij.mll- 
1,9=2,4,...,m (a,y)EQ:; 1,9=2,4,...,.m 
So 


: OPK (Ei, 9) OP k(x,n;) 
OH. HiGe con { ane |v(x)| + er 


ju(y)| 
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6 faeok 
Onno 
We know, 
ju(x)| < ee Va © [xj-2,2;), 
and : 
july) < vy € fay-2.t). 
Therefore, it is not difficult to verify that 
em S Sal gell+ Palgcll+ gligaaesll =m, (3) 


SO 
en = Oh), 


This completes the proof. 


Let fim(£), Gm(X), Kim (x,y) and Uim(a, f(x)), for i = 1,2, are the MHF's 
expansions of f(z), g(x), k;(x,y) and U;(x, f(a)), respectively. According to 
Theorems 1, 2 and expression (1), we have 


Sm(x) ay O(h*) = Im(2) + O(h?) 


4 ([ (kijm(2,y) of O(h*)) (Uijm(y, fim(y) + O(h?)) 4 O(h’)) iy) 


° ( / * (kale, y) + O(K9)) (Wan (us f(y) + O(R3)) + O(83)) 7) 
0 


where x € D. By ignoring the terms included O(h?), we have 


Jt) = amb) + ( J haynes 9) ls fal) 


x ( i * kasmtt, yO am Uh fn(w))a) , (14) 


where x € D. Now, assume the following hypotheses: 
(M1) Suppose that the error of MHFs is denoted by 


Em = ||f — fmll- 


(M2) || fll < M. 
(M3) The nonlinear term U;(z,y) and U2(x,y) satisfies in the Lipschitz 
and linear growth condition such that 


|U1 (a, yr) — Ur (a, yo)| + |U2(a, yr) — U2(a, y2)| < Lys — yo! , 
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where (2, y1), (t,y2) € Dx R, and 


ile, y+ V2le,w)l<L+ly), (aye DxR. 


(M4) Let 
|ki|| <M, 
and 
I|kal| < Mo. 
(M5) Let 
€1m = ||ki — kiyml| < Ci Rh, 
and 


€2,m = ||ko — kaml| < Cgh3, 


where C') and C4 are constants that can be defined as coefficient C in (13). 
(M6) Let 2L(C4h3 + M2)(Mi + C{h3)(L(1 + M) + Cyh3) <1, 
where 
||U1 — Vim] < Crh®, 


and C; is constant that can be defined as coefficient C in (12). 


Theorem 3. Suppose f(x) and fm(a) be the exact and approximate solutions 
of (1) respectively. Also, above assumptions of (M1)-(M6) are satisfied, then 
we have 

Bg = Oh); 


Proof. Assume that w;,m(x) and w;(x) be the approximate and exact solution 
of (8). we define 


Wim(r) = Uim(2, fm(2)), W2,m(x) = U2m(x, fm(2)) ; 


and 
W1m(£) = UilGy Fl) Wm (x) = U2(x, fm(x)) ; 


According to (1) and (14), we have 
f(x) — f(x) = 9(&) — Gm(x) 
+([tulenurtusday) (f° kale udwata, Fun)av 


7 ( | : kim (2, y)Wim(Y, fn(w))dv) ( i : ko,m(2,y)Wam(y, fn(u))av) 


+ O(h*), 


Therefore, 
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«| ( [tales aunts scudiay) — (f° Rarm(tsu)we ls Jv) 


. Cf a,m (2, y)W2m(Ys fn(u))d) 


«| (f rulesaertu seuday) = (f° bayn(ets aden ml Jal) 


+ O(h?). 


So 
Em < em + ||x\|? (Kall |lwil] ||k2w2 — kom w2,mll 


+ lll]? [[R2,mll [v2 nll [itr — Kim 1ml| + O(h*). 


Since x € D, then ||z||? < 1. So 


Em < em + ||K1|| ||wil| ||/kawe — kom we ml 


+ |lko ml [le2,mll R121 — K1,m@1mll + O(h>). (15) 
We have 


where C; and C2 are constants that can be defined as coefficient C in (12). 
Also, we have 


[| Wi mll < [ws — Wim] + [wil] < LEm +CihF+LA+M), (17) 


and 
I[k2,mll < [|k2 — k2,mll + [Rall < Cgh? + Mp . (18) 


Now, according to Theorem 2 and inequalities (14), (15) and assumptions 
(M4)-(M5), we have 


kiwi — Kim Will < |[Kel] we — wWiyml| + llwimll ki — Ki,ml 


< Mi(LEmtCjh?)+Cih3 (LEn+Cih?+L(1+M)). (19) 


From Theorem 1 and inequalities (15)-(17) and assumptions (M3)-(M4), 
we can rewrite (15), as follows 


Em < Ch? + M,L(1+ M) (Lie + CR?) Em + Cah? (Mz + Ch?) 


+LC5h3(1 4 M)) + (Chh° + M2) (LEm + Cih? + L(1-+M)) 
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x (Zan HOPE tC iG PON) AEC a+ M)) + O(n), 


where C is defined in (12). Without loss of generality, ignoring the term 
included E?, and h®, we have 


(C+ ((MLEC) + MpLC})(1-+ M) + MyMa(C1 + C2)) L(+ M))h° 


< 
Em S 1 — 2L(M; + Cih3)(Ma + Ch3)(L( + M) + Che) 


+ O(h?). 
This completes the proof. 


5 Numerical examples 


To illustrate the accuracy and efficiency of proposed method, some ex- 
amples are provided. The algorithms associated with the numerical method 
were performed using Matlab. We have checked that when more points are 
used the accuracy improves significantly. 


Example 1. Consider the following nonlinear QIE [33] 


se)=(@ + 5) 4 (furore) (P Erura). 2 € 0.1, 


with the exact solution f(x) = 2?. 


Table 1 and Figure 1 illustrate the error results for this example. Also, 
we compare the maximum absolute error computed by the present method, 
repeated trapezoidal (RT) method [33] and Adomian decomposition (AD) 
method [33] in Table 2. 


Example 2. Consider the following nonlinear QIE [33] 


payee (FT )+( fA Pena) (fay), en 


where x € [0,1] with the exact solution f(x) = 2. 

Table 3 and Figure 2 illustrate the error results for this example. Also, 
we compare the maximum absolute error computed by the present method, 
RT method [33] and AD method [33] in Table 4. 


Example 3. Consider the following nonlinear QIE 
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Table 1: Absolute error for m = 8, 16,32 of f(a) of Equation (18) 


Nodes x Present method 
m=8 m=16 m=32 

x = 0.0 0 0 0 
x=0.1 8.8569950e-13 = 1.9081958e-17 0 

x = 0.2 2.9361930e-13 1.0984269e-14 6.3490879e-15 
x= 0.3 2.7433902e-11  4.7872678e-12 3.7166104e-13 
x= 0.4 1.4692486e-10  2.1893068e-10 1.1033702e-11 
x = 0.5 1.4339294e-09  1.2865681le-10 1.0435374e-11 
x = 0.6 1.4664141e-07  1.0643045e-08  6.9546968e-10 
x = 0.7 4.2146640e-07  2.6710975e-08  1.5046137e-08 
x = 0.8 6.8445266e-06  3.6807176e-07 5.1551159e-08 
x = 0.9 3.1263755e-06  2.7156875e-06 1.6706609e-07 
x= 1.0 4.2731139e-06  6.9580496e-07  1.1520629¢e-07 


Table 2: Table 2: Comparison of the absolute errors of Example 1 


Methods Maximum error 
RT Method 
m= 10 6.38458E-5 
m = 100 6.30669E-7 
m = 1000 6.30590E-9 
AD Method 
q=5 3.62460E-5 
q = 10 9.23545E-7 
q=15 2.35318E-8 
Present method 
m= 10 1.15435E-5 
m = 100 5.72511E-9 
m = 1000 1.05104E-11 
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Table 3: Absolute error for m = 8, 16,32 of f(a) of Equation (19) 


Nodes x Present method 
m=8 m=16 m=32 

x = 0.0 0) 0 0 
x=0.1 1.7137234e-6 3.3597747e-7 2.0934872c-8 
x = 0.2 6.0344397e-6 3.728595le-7 —1.1336895e-7 
x = 0.3 1.7380694e-5 1.0109144e-6 2.0167842e-7 
x= 04 8.0824578¢e-6 3.4573577e-6  2.0948964e-7 
x = 0.5 1.1099971e-6 7.1694685e-8 5.0373049e-9 
x = 0.6 4.5992070e-5 7.5092211le-6 4.7997155e-7 
x= 0.7 8.4403496e-5 5.4418147e-6 1.7002589e-6 
x = 0.8 1.9036412e-4 1.1238576e-5 2.3184972e-6 
x = 0.9 7.8122715e-5 3.0755442e-5  1.8469227e-6 
x= 1.0 1.9011556e-5 2.0920818e-6 2.8234706e-6 


Table 4: Comparison of the absolute errors of Example 2 


Methods Maximum error 
RT Method 
m= 10 1.07275E-3 
m = 100 8.44338E-7 
m = 1000 8.44337E-9 
AD Method 
q=5 8.44492E-5 
q = 10 8.44338E-7 
q=15 8.44337E-8 
Present method 
m= 10 1.25539E-4 
m = 100 1.27663E-8 
m = 1000 2.35536E-11 
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Absolute error ( logarithmic scale ) 


Absolute error ( logarithmic scale ) 


Figure 2: Absolute errors (on logarithmic scale) for Example 2, with m = 8, 16, 32 


f(e) = 90) + ( [w+ 1) Fay) ( a cos(y)e! dy) , 2€ (0,1, 
(22) 


xz a sin(2x)(1+2z27) xcos(2z) 
1 sin(x) 
( 62 8 4 ) ( ) 


and the exact solution is f(x) = sin(x). 

Table 5 and Figure 3 illustrate the error results for this example. Also, 
we compare the maximum absolute error computed by the present method, 
block-pulse functions (BPFs) method [30] and hat functions (HFs) method 
[28] in Table 6. 
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Table 5: Absolute error for m = 8, 16,32 of f(a) of Equation (20) 


Nodes x Present method 
m=8 m=16 m=32 

x = 0.0 0 0 0 
x=0.1 1.697681 1e-5 3.1978554e-6 1.9245960e-7 
x = 0.2 5.3369972¢e-5 3.0671554e-6 1.0777490c-6 
x = 0.3 1.4115880¢e-4 9.0002882e-6 1.5689730e-6 
x=0.4 5.0416870e-5 2.5604305e-5  1.6783596e-6 
x =0.5 1.5197939e-5 2.3033416e-6  2.9173156e-7 
x = 0.6 2.5006880e-4 3.6295409e-5  2.0906182e-6 
x = 0.7 3.2816702e-4 1.8134154e-5 6.9241732e-6 
x= 0.8 5.4658908¢e-4 3.4864263e-5 6.6704623¢e-6 
x = 0.9 7.9549634e-5 5.0605385e-5  2.3233806e-6 
x= 1.0 1.9488394e-4 4.0158163e-5 6.8004887e-6 


Table 6: Comparison of the absolute errors of Example 3 


Methods Maximum error 
BPFs Method 
m= 8 1.15609E-1 
m= 16 6.36185E-2 
m = 32 3.32452E-2 
HFs Method 
m= 8 2.22432E-2 
m= 16 5.79464E-3 
m = 32 1.51873E-3 
Present method 
m= 8 5.57591 E-4 
m= 16 7.12826E-5 
m = 32 1.10246E-5 


6 Conclusion 


The MHFs method have been proposed for solving nonlinear quadratic inte- 
gral equations. One of the advantages of this method is that the numerical 
solution of the nonlinear QIEs can be converted into a system of algebraic 
equations using the operational matrices. Furthermore, it is proved that 
MHF's method is convergence and the order of convergence is O(h?). The 
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Absolute error ( logarithmic scale ) 


Figure 3: Absolute errors (on logarithmic scale) for Example 3, with m = 8, 16, 32 


proposed method does not need any integration for obtaining the constant 
coefficients hence, it can be applied in a simple and fast technique. The com- 
parison of the obtained results with those based on other methods shows that 
the present method is a powerful mathematical tool for finding the numerical 
solutions of such equations. 
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